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Tunneling delay times of wavepackets in quantum mechanical penetration of 
rectangular barriers have long been known to show a perplexing independence with 
respect to the width of the barrier. This also has relevence to the transmission 
of evanescent waves in optics. Some authors have claimed that in the presence 
of absorption or inelastic channels (which they model by taking a complex barrier 
potential) this effect no longer exists, in that the time delay becomes proportional 
to the barrier width. Taking the point of view that complex potentials imply non- 
Hermitian Hamiltonians and are as such fraught with conceptual pit-falls partic- 
ularly in connection to problems involving time evolution, we have constructed a 
two-channel model which does not suffer from such maladies in order to examine 
this issue. We find that the conclusions arrived at by the earlier authors need to be 
qualified. 



The question of tunneling time [or, 'time spent' by a particle of mass m incident with 
energy(.E) less than the height (Vo) of a barrier of width L in the classically forbidden 
region (E < Vo)] and the analogous situation of evanescent waves in optics, have in recent 
years attracted considerable attention even though interest in such matters go back to 
the forties and fifties of the last century as witnessed by Eisenbud's Ph.D dissertation 1 
with E.P.Wigner at Princeton in 1948 and Wigner's own paper 2 in 1955. David Bohm's 
book on Quantum Theory 3 published in 1951 also throws light on such issues. 

In the present study we shall focus our efforts on one particular aspect of tunnel- 
ing phenomenon [now-a-days referred to as the Hartman 4 -Fletcher 5 effect], namely, the 
counter intuitive conclusion that for a sufficiently opaque barrier tunneling delays are in- 
dependent of the barrier width. By the opaque condition, we imply that the mean energy 
of the incident wave packet is much smaller than the height of the potential barrier. The 
tunneling time (r) is defined in this context by the time taken by the peak of the incident 
wave-packet to traverse the classically forbidden region and to emerge as the transmitted 
wave-packet. Consider the evolution of an incident localized wave-packet described by 
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where Gk {k) is a normalised function of k (say a Gaussian) peaked about the mean 
momentum hko- If this were allowed to evolve freely the peak would travel with a velocity 
hko/m. However with a barrier present an incident plane wave 



\ ^ikx—iEt/h 
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(corresponding to a particle of energy E = where k is the wave vector) would after 
transmission through it emerge as 

-^a T (k)e^ kx - Et ^ h 
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where ar = |a;r|e i<5 is the transmission amplitude (with 5 its associated phase). Accord- 
ingly the transmitted wave-packet would be 

G ko (k)\a T \e ikx - iEt/h+5{k) dk 

Suppose for concreteness, we consider a barrier between x = and x = L, then the 
time(r) at which the peak of the packet will emerge from the barrier at x = L will be 
given using the method of stationary phase by -^[kL — Er/h + 6(h)] = or 

_ fe d8 L 
T ~ n lE + Jhk/m)- [) 

All quantities are to be evaluated at the peak momentum hk and the corresponding 
energy h 2 k 2 /2m. For notational brevity, we suppress the subscript 0. In our discussion, 
this will be taken to be a measure of the tunneling time. 
Specialising to a rectangular barrier, 

V(x) = V for < x < L, (2a) 

V(x) = elsewhere, (26) 

we are to first search for the solution of the Schrodinger equation corresponding to sta- 
tionary states of energy E, 



h 2 d 2 



2m dx 2 



+ V(x) 



ip(x) = Etp(x) (3) 



with the required condition describing reflection and transmission, we have on the two 
sides of the barrier 

ip(x) = e ikx + a R (k)e~ ihx for x < (4a) 

i/>(x) = a T {k)e lkx for x > L, (46) 

where a#(/c) and ar(/c) are the reflection and transmission amplitudes whose square mod- 
ulii give the respective coefficients. The solution within the barrier is given by, 

^(x) = Ae qx + Be~ qx (4c) 
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with q 2 = j!t(V — E). We may as usual proceed to determine A, B, a R and a T by 
matching ip{x) and at the boundaries (x — and x — L) to obtain 

^-(l + |W(- g L) 
(1 + <?) s inh(qL) 

(l-%)exp(qL) , . 

B = ^ FW ; (56) 

(1 + y)sinh{qV) 



and 



4i-exp(— i/cL) 

° T exp(-gL)(l + i-) 2 - exp(qL)(l - i-) 2 6 



and accordingly we arrive at the phase-shift 5. 
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- kL, (7) 



and using eq.(l), we have 
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where the cancellation (in the expression for the time delay) of the contribution arising 
from the term — kL in 5 physically represents the time the wave-packet would have taken 
to traverse the distance from x = to x = L if the potential were absent. It is also 
understood that we must evaluate the delay at the energy Eq. We suppress as mentioned 
earlier the subscript zero for notational convenience. The tunnel time r as given by Eq.(8) 
is plotted as a function of the barrier width L for two values of ^ and shown in Fig.l. As 
must be the case, r — > as L — > 0. But what is more significant is that as L — > oo note 
that tanh(qL) — > 1 and tunnel time r becomes independent of the length L of the barrier. 
This is the counter- intuitive Hartman 4 Effect. The saturated or asymptotic value of r is 
given by 

h m 1 

TaSVmP = - E) = hkq (9) ' 

that is, twice the time taken to traverse the decay distance | in the barrier region if the 



particle were moving with its free velocity — . Explicitly the expression for r is given by 

m 4m m mL 



1 + g 2 \hq 2 h(q 2 — k 2 ) hk 2 Hq 2 sinh(qL)cosh(qL) 
where g = 

It is clearly evident both from the graph and the expression for r that for large L, r 
is independent of L and hence it attains a saturation value. This is what has been 
known as the Hartman 4 -Fletcher 5 effect for decades. What needs to be emphasized is 
that the tunneling time attains the saturation value in two ways depending on the sign of 
k 2 - q 2 = jrr(2E - V ). For k 2 - q 2 > (i.e. E > ^) it monotonically increases to reach 
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Figure 1: The tunnel time (r) is shown as a function of barrier width (L) for two values 
of E/Vq namely 0.1 (dashed curve) and 0.75 (solid curve). 




Figure 2: This graph is for E/V = 0.2. The error in the position of the maxima of the 
hump between the actual value and that calculated from Eq(ll) is about 11%. 



the saturation value and for k 2 — q 2 < (i.e. E < -y) it reaches the saturation value from 
above (there is a hump before it attains saturation). Keeping the next to leading term 
for large L, shows the approach to asymptotia viz. 

(10) 

It is evident from the expression above that there is a exponential approach to the 
asymptotic limit from below if E < ^ (viz. slope is positive with respect to L) and from 
above for E > 

The position of the maximum of the hump (Lq) in the latter case is given by 

1 (hq\ fh\ 10E 2 -9V E + V 2 
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2. Hartman Effect in absorbing media 

For investigating the Hartman 4 -Fletcher 5 effect in presence of an absorbing media, F.Delgado, 
Muga & Rushhaupt 6 and Fabio Raciti & Salesi 7 phenomenologically introduced an imag- 
inary part to the potential in order to describe absorption. This is in the spirit of the 
optical model which is often invoked in nuclear physics when there are inelastic channels. 
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The relative phase-shift between the transmitted and the incident wave-packet in such a 
situation is obtained by just putting the complex potential form in place of Vq in q and 
inserting it in eq(7). If the potential is given by V = Vq(1 — iX) within the range L and 
zero elsewhere, then the transmission coefficient (a^) is given by 

Hence in the large L limit, ax will take the form of 

ttT = e^(l-f) 2 (13) 
For complex potentials and taking only the first order term in A, it is seen that 

iXVo 



q = q °V 2(V -E) 
where q = [^{V — E)} 1 / 2 . Hence, we have 



q 90 V 2(V„ - E) 
which is readily expressed in a polar form as 




2q {V Q - E) 

Hence ax in presence of the complex potential will be 

( 1 + Wo \ -ikL-q (l- 2 ^° E) 
\ L ^ 2(V -E) J e 



M] > ^ 2 e ^- l (w^) (14) 
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ax = 
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(15) 



The only L dependence in the phase-shift comes from the numerator of eqn(15) and the 
dependence is of the form of 

Hence the tunneling time r for large L is given by 

XV k 



2(V -E)q 



L + i (16) 



where ^ is the L-independent factor that comes from the denominator of eqn(15). It is 
evident from the expression of r that the tunneling time for large L is not independent 
of L and hence we no longer have the Hartman-Fletcher effect even for weak absorption 
(small A). Optical experiments by Nimtz, Spieker & Brodowsky 8 confirms the absence of 
Hartman-Fletcher effect only for strong absorption and the absence of this effect is not 
claimed for weak absorption. 
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Furthermore, one may express doubts regarding the conclusion on tunneling times in 
the presence of inelastic channels based on the phenomenological device of introducing 
complex potentials (through the paradigm of the optical model) particularly when we are 
dealing with questions regarding time evolution. Non-unitary time-evolution using a non- 
hermitian Hamiltonian could be questionable. If randomness or disorder in the system 
is invoked then one must use a density matrix approach to the netire problem. We put 
forward a two-channel formulation for this problem and analyze the matter of tunneling 
time in this format where unitarity is properly maintained. 

To obtain a simple realisation of such a situation, consider the elastic scattering of a 
particle by a target system describable in terms of a rectangular potential Vq of width L. 
Suppose that the target has an excited state of energy A above its ground state. When 
the projectile energy E is above the threshold A then, apart from the elastic scattering 
the inelastic process is also possible. One may however describe the situation in terms 
of two channels, the elastic and the inelastic. The latter channel corresponds to the 
scattering of the projectile of energy E > A from the excited state of the system through 
a rectangular potential Vj say of same width L, while the two channels are coupled to each 
other via a potential V c for simplicity, taken, to be that of same width L. With ip(x) and 
<f)(x) denoting the wave-functions in the elastic and the inelastic channels respectively, the 
system is governed by the coupled Schrodinger equations 

ijj(x) + V c <j>{x) = E^j(x) (17a) 

4>{x) + V c ip(x) = (E — A)(f)(x) (17b) 

The wave-functions for the elastic channel and the inelastic channel in the region x < L 
is given as 

rj}{x) = e ikx + a R (k)e~ ikx (18a) 

4>(x) = Re~ ik ' x (18b) 

where k 2 = and k'2 = 2m (^ 2 ~ A ) while an(k) and R are the respective reflection 
amplitudes in the elastic and inelastic channel. Similarly the wave-functions in the region 
x > L is 

i)(x) = a T (k)e lkx (19a) 

<P(x) = Te lk ' x (196) 

where a^ and T are the respective transmission amplitudes in the elastic and inelastic 
channels. Note here that naturally there are two outgoing waves to the left and to the 
right due to the additional inelastic channel as physically required. The wave-functions 



in the two channels in the intermediate region are 

i/>(x) = (Be ax + Ce- ax )sin(6/2) + (Fe Px + Ge- px )cos(9 /2) (20a) 
0(x) = (Be ax + Ce- ax )cos(9/2) - (Fe? x + Ge^ x )sin(9/2) (20b) 

where 

a 2 = V ° + V >~ 2E + A + -Vj- A) 74 + V c 2 (21a) 



2m dx 2 
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P 2 = 



V + Vr - 2E + A 
2 



y/(V -Vj- Ay/4 + V? 



(216) 



= 



(21c) 



Vfc - Vi - A 



The phase difference between the incident and the transmitted waves in the elastic channel 
can be easily computed from ax evaluated by demanding the continuity of the wave- 
function and its derivative at the boundaries and hence the tunneling time in the elastic 
channel is expressed as a function of a and f3. The two channels, being coupled by a 
coupling potential V c , must admit evanescent wave-functions inside the potential barriers 
of both the channels for Hartman-Fletcher effect to be present in the elastic channel. For 
the wave-function to be evanescent inside, a and (3 must be real. Now 



For weak coupling (i.e., low V c ) and 'opaque condition' in the inelastic channel, i.e, Vj > 
E — A, we have a/3 > and real. If the opaque condition in both the channels are to be 
respected along with weak coupling, then both a and (3 cannot simultaneously be purely 
imaginary. Hence both a and (3 are real as a result of which evanescent wave occurs inside 
the barriers in the two-channels. However, if Vr < E — A, then aj3 is imaginary and both 
or either of them must be imaginary resulting in the absence of evanescent wave inside 
the barrier in the elastic channel. In general, Hartman-Fletcher effect is expected if 



Instead, if an attractive potential is considered in the inelastic channel, which is physically 
not forbidden, we have Vr = —V® (say), then 



If the inelastic channel is to be kept open, i.e., E > A, then a(3 is imaginary even for a 
weak coupling. So either or both of a, (3 are imaginary, as a result of which no evanescent 
wave, and hence the Hartman-Fletcher effect is expected. 

We have thus demonstrated that the independence of the tunneling time with respect 
to the width (L) of a barrier for a real potential with large L, claimed to be violated if 
inelastic channels are present, is a statement which needs to be qualified. We show that 
if instead of trying to capture inelasticity through the device of a complex potential, we 
treat the process through a multichannel formalism (more appropriate for the study of 
time evolution as compared to approaches using on-Hermitian Hamiltonians), then the 
outcome depends on some details of the channels and the strength of the coupling between 
them. 

*present address: Max-Plank-Institut fur Physik komplexer Systeme, Nothnitzer Stra/3e 38, D-01187 
Dresden, Germany. 



a(3= y/(V T -E + A)(V -E) -V c 2 . 



(22) 



(Vj - E - A)(V - E) > V c 2 



(23) 




(24) 
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